n=1e4, d=16^7, p = 1 - exp[(-n*(n-1))/(2*d)] , q = 1 - [(d-1)/d]^n

Uh oh! Wolfram|Alpha doesn't run without JavaScript.

Please enable JavaScript. If you don't know how, you can find instructionshere.Once you've done that, refresh this page to start using Wolfram|Alpha.